Forecasting the Cosmological Constraints with Anisotropic Baryon Acoustic 

Oscillations from Multipole Expansion 



O 
(N 

in 

(N 



o 

u 

6 



> 
m 



Atsushi Taruya^'^, Shun Saito'^''', Takahiro Nishimichi^ 
^Research Center for the Early Universe, School of Science, 
The University of Tokyo, Bunkyo-ku, Tokyo 113-0033, Japan 
'^Institute for the Physics and Mathematics of the Universe, 
The University of Tokyo, Kashiwa, Chiba 277-8568, Japan 
^Department of Physics, The University of Tokyo, 113-0033, Japan and 
'''Department of Astronomy, University of California Berkley, California 94720, USA 

(Dated: January 26, 2011) 

Baryon acoustic oscillations (BAOs) imprinted in the galaxy power spectrum can be used as a 
standard ruler to determine angular diameter distance and Hubble parameter at high redshift galax- 
ies. Combining redshift distortion effect which apparently distorts the galaxy clustering pattern, 
we can also constrain the growth rate of large-scale structure formation. Usually, future forecast 
for constraining these parameters from galaxy redshift surveys has been made with a full 2D power 
spectrum characterized as function of wavenumber k and directional cosine fi between line-of-sight 
direction and wave vector, i.e., P{k, Here, we apply the multipole expansion to the full 2D power 
spectrum, and discuss how much cosmological information can be extracted from the lower-multipole 
spectra, taking a proper account of the non-linear effects on gravitational clustering and redshift 
distortion. The Fisher matrix analysis reveals that compared to the analysis with full 2D spectrum, 
a partial information from the monopole and quadrupole spectra generally degrades the constraints 
by a factor of ~ 1.3 for each parameter. The additional information from the hexadecapole spectrum 
helps to improve the constraints, which lead to an almost comparable result expected from the full 
2D spectrum. 



I. INTRODUCTION 

Baryon acoustic oscillations (BAOs) imprinted on the 
clustering of galaxies are now recognized as a powerful 
cosmologicalprobe to trace the expansion history of the 
Universe [ll-Q- particular, the BAO measurement via 
a spectroscopic survey can provide a way to simultane- 
ously determine the angular diameter distance Da and 
Hubble parameter H at given redshift of galaxies through 
the cosmological distortion, known as Alcock-Paczynski 
effect (e.g., [1-Q). Further, measuring the clustering 
anisotropics caused by the redshift distortion due to the 
peculiar velocity of galaxies, we can also probe the growth 
history of structure formation (e.g., |9l-[l2|). characterized 
by the growth-rate parameter / = dhiD/dlna, with 
quantities D and a being linear growth factor and the 
scale factor of the Universe, respectively. 

With the increased number of galaxies and large sur- 
vey volumes, on-going and future spectroscopic galaxy 
surveys such as Baryon Oscillation Spectroscopic Sur- 
vey (BOSS) ^13i], Hobby- Eberly Dark Energy Experi- 
ment (HETDEX) [l^ , Subaru Measurement of Imaging 
and Redshift equipped with Prime Focus Spectrograph 
(SuMIRe-PFS), and EUCLID/JDEM [15, .16,] aim at pre- 
cisely measuring the acoustic scale of BAOs as a stan- 
dard ruler. These surveys will cover the wide redshift 
ranges, 0.3 ^ z < 3.5, and provide a precision data of 
the redshift-space power spectrum with an accuracy of a 
percent level over the scales of BAOs. 

In promoting these gigantic surveys, a crucial task is a 
quantitative forecast for the size of the statistical errors 
on the parameters Da, H and / in order to clarify the 



scientific benefits as well as to explore the optimal survey 
design. The Fisher matrix formalism is a powerful tool to 
investigate these issues, and it enables us to quantify the 
precision and the correlation between multiple parame- 
ters (d, 0, 113, [Hi, especially for measuring Da, H and 
/). So far, most of the works on the parameter forecast 
study have focused on the potential power of the BAO 
measurements, and attempt to clarify the achievable level 
of the precision for the parameter estimation. For this 
purpose, they sometimes assumed a rather optimistic sit- 
uation that a full shape of the redshift-space power spec- 
trum, including the clustering anisotropies due to the 
redshift distortion, is available in both observation and 
theory. 

In this paper, we are particularly concerned with the 
parameter estimation using a partial information of the 
anisotropic BAOs from a practical point-of-view. In red- 
shift space, the power spectrum obtained from the spec- 
troscopic measurement is generally described in the two 
dimension, and is characterized as functions of k and ^, 
where k is the wavenumber and n is the directional co- 
sine between the line-of-sight direction and k [i^ . While 
most of the forecast study is concerned with a full 2D 
power spectrum, the multipole expansion of redshift- 
space power spectrum has been frequently used in the 
data analysis to quantify the clustering anisotropies. De- 
noting the power spectrum by P(fc,/i), we have 
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with the function Vt being the Legendre polynomials. Al- 
though the analysis with full 2D spectrum will definitely 
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play an important role as improving the statistical sig- 
nal, most of the recent cosmological data analysis has fo- 
cused on the angle-averaged power spectrum = 0), i.e., 
monopole spectrum, and a rigorous analysis with full 2D 
spectrum is still heavy task due to the time-consuming 
covariance estimation (e.g., [l^-[2l|). 

In linear theory, the redshift-space power spectrum is 
simply written as P(fc,/i) = (1 -|- /3 /i^)^Pgai(fc), where 
(3 — f/b with b being the linear bias parameter, and 
Pgai is the galaxy power spectrum in real space l22 - l24l |. 
Then, the non- vanishing components arises only from the 
monopole = 0), quadrupole {£ = 2) and hexadecapole 
spectra {£ — 4). That is, cosmological information con- 
tained in the £ = 0, 2 and 4 moments is equivalent to 
the whole information in the full 2D power spectrum. 
Observationally, however, this is only the case when we 
a priori know the cosmological distance to the galaxies. 
The Alcock-Paczynski effect can induce non-trivial clus- 
tering anisotropics, which cannot be fully characterized 
by the lower multipole spectra, in general. Further, in re- 
ality, linear theory description cannot be adequate over 
the scale of the BAOs, and the non-linear effects of the 
redshift distortion as well as the gravitational cluster- 
ing must be accounted for a proper comparison with ob- 
servation. These facts imply that non-vanishing multi- 
pole spectra higher than £ > 4 generically appear, and a 
part of the cosmological information might be leaked into 
those higher multipole moments. An important question 
is how much amount of the cosmological information can 
be robustly extracted from the lower multipole spectra 
instead of the full 2D spectrum. In the light of this, 
Ref. (H recently examined a non-parametric method to 
constrain Da and H from the monopole and quadrupole 
spectra, and numerically estimate the size of errors (see 
also Ref. [1^ for the estimation of growth-rate parame- 
ter). 

Here, as a complementary and comprehensive ap- 
proach, we will investigate this issue based on the Fisher 
matrix formalism, and derive the useful formulae for pa- 
rameter forecast using the multipole power spectra. We 
then explore the potential power of the lower multipole 
spectra on the cosmological constraints, particularly fo- 
cusing on the parameters Da, H and /. To do so, we 
consider the Figure-of-Merit (FoM) and Figure-of-Bias 
(FoB) for these parameters, and investigate their de- 
pendence on the assumptions for the number density of 
galaxies, the amplitude of clustering bias, the maximum 
wavenumber used for the parameter estimation. 

In Sec. nil we present the Fisher matrix formalism for 
cosmological parameter estimation from the multipole 
power spectra. Sec. IIIII deals with the model of redshift- 
space power spectrum and the assumptions used in the 
Fisher matrix analysis. Then, in Sec. IIV| the results for 
FoM and FoB are shown, and the sensitivity of the re- 
sults to the assumptions and choice of the parameters is 
discussed in greater details. Finally, Sec. |V] briefly sum- 
marize our present work. 

Throughout the paper, we assume a flat Lambda cold 



dark matter (CDM) model, and the fiducial model pa- 
rameters are chosen based on the five-year WMAP results 
[2^: = 0.279, 17a = 0.721, f^b = 0.0461, h = 0.701, 
n, = 0.96, A, = 2.19 x 10"^. 



II. FISHER MATRIX FORMALISM 

In this section, we present the basic formulae for Fisher 
matrix analysis in estimating the statistical error and 
systematic biases for cosmological parameters from the 
multipole power spectra. 

Let us first derive the expression for Fisher matrix rel- 
evant for the power spectrum analysis. The definition of 
the Fisher matrix is given by 
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where 9i denotes the parameter, and the quantity C is 
the likelihood function. For the parameter estimation 
study with the multipole spectrum, Pi{k), the likelihood 
function is usually taken in the form as 



C oc exp 



(3) 



where we define 



APt{k) = Peik)-Pi{k), 
C"'{k^,kn) = (AP,(fc„OAPr(fc„)). 

The quantities Pi{k) and Pi{k) respectively denote the 
observed estimate and theoretical template for the mul- 
tipole power spectrum. 

Substituting Eq. Q into the definition the leading- 
order evaluation of the Fisher matrix leads to (e.g., [23, 
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where we have assumed that the covariance is approx- 
imately characterized by the Gaussian statistic, and is 
written as C"'(fcm,fc„) = Gov" (fc„)5„„. 

Adopting the power spectrum estimation by Ref. [2^ , 
the analytic expression for the quantity Gov^^ (fc„) can 
be found in Ref. [s^l [see Eq. (25) of their paper]: 



Jd^r n(r)2[l + n(r)P(s)(fc„,^)] 



(5) 



with Pi{fJ-) being the Legendre polynomiaUSO] . The 
quantity Vn is the volume element of the thin shell in 
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the Fourier space, i.e., Vn = 47r^fc^dfc„/(27r)'^, which cor- 
responds to AT4/(27r)^ in the notation of Ref. (soj . 

Now, to simphfy the formula, we consider the homoge- 
neous galaxy samples, which implies n{r) = n = const. 
In this case, the denominator in the integrand of Eq. ([5]) 
is simplified as 



^V, |P(fc,Ai) + -| , (6) 

where Vs denotes the survey volume. Then, taking the 
continuum limit, the expression for Fisher matrix can be 

recast as 
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with the reduced covariance matrix Gov (fc) given by 
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Here, the range of integration [k„ 
sen through the survey properties and/or limitation of 
theoretical template, and, in particular, the minimum 

1 /3 

wave number is limited to 27r/Vi . 

Eq. ([7]) with ([S]) is the formula for the Fisher matrix 
used in the parameter estimation with multipole power 
spectra. This can be compared with the standard for- 
mula for fuU 2D power spectrum (e.g., [1,0, [28,]): 
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That is, the full 2D information obtained through the in- 
tegral over directional cosine [i, in Eq. ([9|) is replaced with 
the summation over all multipoles in the new formula 
([7]). Thus, truncating the summation at a lower multi- 
pole generally leads to the reduction of the amplitude in 
Fisher matrix, and as a result, the statistical errors of the 
parameter Q j margin alized over other parameters, given 
by A6'i = ^ {i^^ijii, is expected to become larger. 

The Fisher matrix formalism also provides a simple 
way to estimate the biases in the best-fit parameters 
caused by the incorrect template for the multipole power 
spectra P™'^^{k). To derive the formula for system- 
atic bias, we replace the template power spectrum Pe{k) 
in the likelihood function ([3]) with the incorrect one 
pwiong^^-j^ We denote this likelihood function by As- 
suming that the size of the biases are basically small, the 



(biased) best-fit values can be estimated from the ex- 
tremum of the Likelihood function C by expanding the 
expression of the extremum around the fiducial parame- 
ters: 
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where the quantities with subscript stand for the one 
evaluated at the fiducial parameters, and the 59i means 
the deviation of the best-fit value from the fiducial pa- 
rameter. Then, taking the ensemble average of the above 
expressions and using the definition of the Fisher matrix, 
we obtain 



(11) 



where the Fisher matrix F^^ is the same one as given by 
Eq. d?]), but is evaluated using incorrect power spectra 
pwrong(^^_ The vector is 
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Here, the multipole power spectrum P^^^{k) denotes 
the systematic difference between correct and incorrect 
model of multipole power spectra, P^^^{k) — P™'^^{k) — 
P^'^^°{k). In deriving the above expression, we have used 
the fact that the extremum of the likelihood function 
vanishes only when the correct template for the multi- 
pole power spectrum is applied. 

Notice that similar but essentially different formula for 
systematic biases is obtained in the cases using the full 
2D power spectrum. It is formally expressed as Eq. pT|) . 
but the Fisher matrix F/^ is now replaced with Eq. ^ 
evaluated using the incorrect 2D spectrum P^''°"s(/j^ ^). 
Further 

for the full 2D spectrum (e.g., [3ll |3 



the vector Sj should be replaced with the one 
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Finally, all the formulae derived in this section ignore 
the non-Gaussian contributions to the likelihood and co- 
variances, which would be sometimes important in prac- 
tice. The extension of the formulae to include the non- 
Gaussian contributions is straightforward, and will be 
considered elsewhere. For the effects of non-Gaussian 
contributions to the parameter estimation study espe- 
cially focusing on BAOs, several works have been recently 
done based on the numerical and analytical treatments 

[Hill HI. 
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FIG. 1: Monopole (top) and quadrupole (bottom) moments of 
matter power spectra in redshift space at z = 1. The results 
are divided by the smooth reference spectrum, no-wiggid 
and are compared with the N-body results (symbols) taken 
from the wmap5 simulations of Ref. 35]. The reference spec- 
trum no-wiggic calculated from the no- wiggle approxima- 
tion of the linear transfer function [36l ] with the linear theory 
of the Kaiser effect taken into account. Solid and dot-dashed 
lines represent the results of improved PT calculations based 
on the model of redshift distortion (|14p , but the terms A and 
B are ignored in the dot-dashed lines. In both cases, the 
one-dimensional velocity dispersion CTv was determined by fit- 
ting the predictions to the N-body simulations, using the data 
below the wavenumber indicated by the vertical arrow. The 
best-fit values of ctv are CTv — 395 kms~^ and 285 kms~^, with 
and without the A and B terms, respectively. 



III. MODEL AND ASSUMPTIONS 

Given the formulae for Fisher matrix analysis, we now 
move to the discussion on the parameter forecast study 
using the multipole power spectra, and compare the re- 
sults with those obtained from the full 2D spectrum. Be- 
fore doing this, in this section, we briefly describe the 
model and assumptions for redshift-space power spec- 
trum relevant for spectroscopic measurement of BAOs. 

In redshift space, clustering statistics generally suffer 
from the two competitive effects, i.e., enhancement and 
suppression of clustering amplitude, referred to as the 
Kaiser and Finger-of-God effects, respectively. While the 
Kaiser effect comes from the coherent motion of the mass 
(or galaxy), the Finger-of-God effect is mainly attributed 
to the virialized random motion of the mass residing at 
a halo. On weakly non-linear regime, a tight correlation 
between velocity and density fields still remains, and a 
mixture of Kaiser and Finger-of-God effects is expected 
to be significant. Thus, a careful treatment is needed for 



accurately modeling anisotropic power spectrum. 

Recently, we have presented an improved prescription 
for matter power spectrum in redshift space taking ac- 
count of both the non-linear clustering and redshift dis- 
tortion [32[. Based on the perturbation theory calcu- 
lation, the model can give an excellent agreement with 
results of N-body simulations, and a percent level preci- 
sion is almost achieved over the scales of our interest on 
BAOs. The full 2D power spectrum of this model is very 
similar to the one proposed by Ref. [s^l , but includes the 
corrections: 

+ f / Peeik) + A{k, fi; /) + B{k, ^; /) } (14) 

with the quantity / being the growth-rate parameter. 
Here, the power spectra Pss, Pee and Pse denote the auto 
power spectra of density and velocity divergence, and 
their cross power spectrum, respectively. The velocity di- 
vergence 9 is defined hy 9 = —Vv/{aHf). The qua ntity 
(Tv denotes the one-dimensional velocity dispersion [5 ij. 
and the exponential prefactor characterizes the damping 
behavior by the Finger-of-God effect. For the purpose to 
model the shape and structure of BAOs in power spec- 
trum, (Tv may be treated as a free parameter, and deter- 
mine it by fitting the predictions to the observations. 

A salient property of the model is the presence 
of the terms A and B, which represent the higher-order 
couplings between velocity and density fields, usually ne- 
glected in the phenomenological models of redshift dis- 
tortion. The explicit expressions for these terms are de- 
rived based on the standard treatment of perturbation 
theory, and the results are presented in Ref. [s^l. A de- 
tailed investigation in our previous paper [32*1 reveals that 
the corrections A and B can give an important contribu- 
tion to the acoustic structure of BAOs over the scales 
k ^ 0.2hMpc~^ , which give rise to a slight uplift in the 
amplitude of monopole and quadrupole spectra. With 
the improved treatment of the perturbation theory to 
compute Pss, Pee and Pse (e.g., [H,!!!]), the model p4)) 
can give a better prediction than the existing models of 
redshift distortion. Fig. [T] plots the illustrated example 
showing that the model (|14l) reproduces the N-body re- 
sults of monopole and quadrupole spectra quite well, and 
the precision of the agreement between prediction and 
simulation reaches a percent-level. Hence, in this paper, 
we adopt the model as a fiducial model for matter 
power spectrum in redshift space. 

Note that the model generically produces the non- 
vanishing higher multipole spectra oi £ > 4, due to the 
damping factor, e~'^^^^'^^^ . Furthermore, the correc- 
tions A and B are expanded as power series of /i, which 
include the powers up to /i^ for the A term, /x* for the 
B term. This indicates that the corrections addition- 
ally contribute to the higher multipoles, at least, up to 
£ = 8. In this sense, the model provides an inter- 
esting testing ground to estimate the extent to which the 
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useful cosmological information can be obtained from the 
lower-multipole spectra. 

Then, assuming the hnear galaxy bias in real space, 
Sgai = Wmass, the redshift-space power spectrum for 
galaxies becomes 

Pgai ik,^l)= e-(^^ f^-^" [Pss (fc) + 2 /? ^? Pse (k) 

+ 13^ / Pgeik) + bA{k, /i; /3) + B{k, ^i; /?)} (15) 

with /? — f/b. The linear deterministic bias may 
be too simplistic assumption, and the effects of non- 
linearity and stochasticity in the galaxy bias might be 
non-negligible js^-Hlj. Our primary concern here is the 
qualitative aspects of the parameter estimation using the 
multipole spectra, based on a physically plausible model 
of redshift distortion. Since the galaxy bias itself does 
not produce additional clustering anisotropics, we sim- 
ply adopt the linear bias relation for illustrative purpose. 

Finally, notice that in addition to the clustering 
anisotropics caused by the peculiar velocity of galaxies, 
the observed galaxy power spectrum defined in comov- 
ing space further exhibits anisotropics induced by the 
Alcock-Paczynski effect. This is modeled as 

where the quantity Pgai(9, i^) at the right- hand-side repre- 
sents the template for the redshift-space power spectrum 
in the absence of cosmological distortion, i.e., Eg. dTSI) . 
The comoving wavenumber k and the directional cosine 
fj, measured with the underlying cosmological model are 
related to the true ones q and ly by the Alcock-Paczynski 
effect through (e.g., 0,141, El) 




The quantities DA,fid and iJfid are the fiducial values of 
the angular diameter distance and Hubble parameter at 
a given redshift slice. 



IV. RESULTS 

In what follows, for illustrative purpose, we consider 
the hypothetical galaxy survey of the volume Vs — 
4/i~^Gpc^ at z = 1, and examine how well we can con- 
strain the distance information and growth-rate parame- 
ter. Da, H, and /, from the low- multipole power spectra. 



We set the number density of galaxies, linear bias param- 
eter and velocity dispersion to n = 5 x 10~'*/i'^Mpc~^, 
b — 2 and CTv = 395km s~^. These values are used in 
the Fisher analysis as a canonical setup, but we also ex- 
amine the variants of these parameter set to study the 
sensitivity of the forecast results. Note that the depth 
and the volume of the survey considered here roughly 
match those of a stage III class survey defined by the 
Dark Energy Task Force (DETF) [H]. 

To compute the Fisher matrix adopting the model of 
redshift-space power spectrum, Eq. (|15|) . we just follow 
the procedure in Ref. [33 to calculate the redshift-space 
power spectra. That is, we use the improved PT devel- 
oped by Ref. (ssl . [isf to account for a dominant contribu- 
tion of the non- linear gravity to the power spectra Pss, 
Pse and Pge, and to adopt standard PT for small but non- 
negligible corrections of A and B terms. Detailed com- 
parison with N-body simulations f32','35'l showed that this 
treatment can work well, and in our fiducial set of cos- 
mological parameters, the model can give a percent-level 
precision at least up to the wavenumber k < 0.2/iMpc~"'^ 
at z = 1. 

Number of free parameters in the subsequent Fisher 
analysis is five in total, i.e.. Da, H, and /, in addition to 
the parameters b and a^- Other cosmological parameters 
such as rim or fib are kept fixed. We assume that the 
cosmological model dependence of the power spectrum 
shape is perfectly known a priori from the precision CMB 
measurement by PLANCK [4^. The influence of the 
uncertainty in the power spectrum shape is discussed in 
Sec. live 21 in detail. 



A. Two-dimensional errors 

As a pedagogical example, let us first examine how 
the lower-multipole spectra can constrain the parameters 
Da, H, and /. Fig. [5] shows the two-dimensional contour 
of the 1-cr (68% C.L.) errors on {Da,H) (bottom-left), 
{Da,/) (top-left), and (/, iJ)-planes (bottom right). 
Here, the Fisher matrix is computed adopting the model 
of redshift-space power spectrum (jl5l) up to fcmax = 
0.2/iMpc-i. 

The magenta solid and cyan dashed lines respectively 
represent the constraints coming from the monopole (Pq) 
and quadrupole (P2) power spectrum alone. As antici- 
pated, only the single multipole spectrum cannot pro- 
vide useful information to simultaneously constrain Da, 
H , and /. In particular, for the constraints on Da and 
H, there appear strong degeneracies, and the error el- 
lipses are much elongated and inclined. These behaviors 
are basically deduced from the Alcock & Paczynski ef- 
fect, and are consistent with the facts that the monopole 
spectrum is rather sensitive to the combination {D\/H), 
while the quadrupole spectrum is sensitive to {Da H) 
(e.g., d). On the other hand, combining monopole 
and quadrupole greatly improves the constraints (indi- 
cated by blue, outer shaded region) not only on Da and 
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iJ, but also on growth-rate parameter /. This is be- 
cause the degeneracies between the parameters Da and 
H constrained by the monopole differ from that by the 
quadrupole, and thus the combination of these two spec- 
tra leads to a substantial reduction of the size of error 
ellipses. Further, the growth-rate parameter is propor- 
tional to the strength of redshift distortion, and can be 
determined by the quadrupole-to-nionopole ratio. Al- 
though the measurement of the galaxy power spectrum 
alone merely gives a constraint on j3 = f /h, provided 
the accurate CMB measurement for power spectrum nor- 
malization, we can separately determine the growth-rate 
parameter. Note that the combination of monopole and 
hexadecapole spectra also provides a way to determine 
the growth-rate parameter (red shaded region), although 
the error on / is a bit larger due to the small amplitude 
of hexadecapole spectrum. 

For comparison. Fig. [5] shows the forecast constraints 
obtained from the full 2D power spectrum (green, in- 
ner shaded region). Further, we plot the results com- 
bining the monopole and quadrupole spectra, but ne- 
glecting the covariance between £ = and £ = 2, i.e., 

02 20 

Gov — Gov = (blue, dotted lines). Glearly, us- 
ing a full 2D shape of the redshift-space power spec- 
trum leads to a tighter constraint, and the area of the 
two-dimensional error is reduced by a factor of 1.6 — 18, 
compared with the constraints from the monopole and 
quadrupole spectra. The results indicate that the con- 
tribution of the higher multipoles is very important, and 
the additional information from quadrupole and hexade- 
capole spectra, each of which puts a different parameter 
degeneracy, seems to play a dominant role in improving 
the constraints. On the other hand, for joint constraints 
from the monopole and quadrupole, a role of the covari- 

^02 20 

ance Gov or Gov seems less important, and one may 
naively treat monopole and quadrupole power spectra as 
statistically independent quantities. However, these re- 
sults are partially due to the properties of the galaxy 
samples characterized by several parameters, and may 
be altered in different assumptions or survey setup. This 
point will be investigated in some details in next subsec- 
tion. 



B. Figure-of- Merit 

We here study the dependence of galaxy samples or 
survey setup on the forecast results for parameter con- 
straints. To do this, it is useful to define the Figure-of- 
Merit (FoM): 

FoM = ^ (19) 
V dctF~^ 

— 1 

where the matrix F is the 3x3 sub-matrix, whose 
elements are taken from the inverse Fisher matrix 
associated with the parameters Da, H, and /. 




FIG. 2: Two dimensional contours of 1-cr (68%CL) er- 
rors on {Da,H) (bottom-left), {Da, f) (top-left), and {f,H) 
(bottom-right), assuming a stage-Ill class survey with Vs = 
4/i~^Gpc^ at z = 1. In each panel, magenta solid and cyan 
dashed lines respectively indicate the forecast constraints 
coming from the monopole (Po) and quadrupole (P2) spec- 
trum alone, while the blue and red shaded region represent 
the combined constraints from Po and P2, and Po and P4, 
respectively. The inner green shaded region is the results 
coming from the full 2D spectrum. As a reference, blue dot- 
ted contours show the results combining both Po and P2, but 

(incorrectly) neglecting the covariance between monopole and 

02 ^20 

quadrupole spectra, i.e.. Gov = Gov — 0. 



The FoM quantifies the improvement of the parame- 
ter constraints, and inversely proportional to the prod- 
uct of one-dimensional marginalized errors, i.e., FoM cx 
l/{a{DAMH)aif)}. 

Fig. [3] shows the dependence of FoM on the properties 
of the galaxy samples characterized by the number den- 
sity Ug (top-right), bias parameter b (bottom- left), and 
one-dimensional velocity dispersion (bottom- right). 
Also, in top-left panel, we show the FoM as a function of 
maximum wavenumber fcmax used in the parameter esti- 
mation study. Note that in plotting the results, the other 
parameters are kept fixed to the canonical values. The 
upper part of each panel plots the three different lines, 
and shows how the FoM changes depending on the choice 
or combination of power spectra used in the analysis: 
combining monopole (Pq) and quadrupole (P2) spectra 
(magenta, dot-dashed); combining three multipole spec- 
tra, Po, P2 and P4 (blue, long-dashed); using full 2D 
spectrum P{k,ix) (black, solid). On the other hand, the 
lower part of each panel plot the ratio of FoM normalized 
by the one for the full 2D spectrum. 
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FIG. 3: Figure-of- merit (FoM) on the parameters 7?a, H, and / defined by Eq. (|19|l . as functions of fcmax (top-left), rigai (top- 
right), 6 (bottom left), and ctv (bottom right), assuming a hypothetical gala:xy survey at z = 1 with volume Va = 4h~^Gpc^. 
In each panel, solid lines are the results obtained from the full 2D power spectrum, while the dashed and dot-dashed lines 
represent the FoM from the combination of the multipole spectra (dot-dashed: Po & P2, dashed: Po, Po, & Pa)- The bottom 
panels show the ratio of FoM normalized by the one obtained from the full 2D spectrum. Note that except the parameter along 



the horizontal axis, the fiducial values of the model parameters are set to fcn 
and (Tv = 3.95 /i~^Mpc, indicated by the vertical dotted lines. 



0.2/iMpc-S Ug = 5 X 10-*/^^Mpc-^ b = 2, 



In principle, using the full 2D spectrum gives the tight- 
est constraints on Da, H , and /, but an interesting point 
here is that almost equivalent FoM to the one for the 
full 2D spectrum is obtained even from a partial infor- 
mation with the lower- multipole spectra Pq, P2 and P4. 
This is irrespective of the choice of the parameters for 
galaxy samples. Although the result may rely on the 
model of redshift distortion adopted in this paper, re- 
calling the fact that the non-vanishing multipole spec- 
tra higher than £ ^ 6 arise only from the non-linear 
effects through the gravitational evolution and redshift 



distortion, the cosmological model dependence encoded 
in these higher multipoles is expected to be very weak, 
partly due to the low signal-to-noise ratio. In this sense, 
the result in Fig. [3] seems reasonable. 

Now, turn to focus on the FoM from the combination 
of Po and P2. Fig. [3] indicates that except for the case 
varying the bias 6, the resultant FoM shows a mono- 
tonic dependence on the parameters. As a result, the 
ratio of FoM shown in the lower part of the panels is 
nearly constant around 0.4 — 0.6. As for the variation 
of bias parameter, the non-monotonic dependence of the 
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FIG. 4: Correlation coefficient for tlie covariance, rcov = 

COT°'V[Cw°'°COT^'Y^^ as function of /3 = //b. Tlie plot- 
ted results are obtained based on tlie linear theory, in which 
the coefficient rcov depends on the power spectrum ampli- 
tude relative to the shot-noise contribution, rig P, as well as 
fi. The solid, long-dashed, short-dashed, and dot-dashed lines 
respectively indicate the results with rig P = 1, 2, 5, and 10. 



FoM is basically explained by the two competitive effects. 
That is, as increasing 6, while the power spectrum ampli- 
tude increases and signal-to-noise ratio is enhanced, the 
clustering anisotropics due to the redshift distortion con- 
trolled by the quantity /? are gradually reduced. Hence, 
for some values of 6, FoM becomes maximum. A notice- 
able point is that the ratio of FoM for the monopole and 
quadrupole becomes gradually increased as the clustering 
bias becomes large. At 5 4, the ratio of FoM reaches at 
0.8, indicating most of the cosmological information con- 
tained in the hexadecapole and higher multipoles is lost, 
and signals coming from the monopole and quadrupole 
spectra becomes dominated. 

The reason for this behavior is presumably due to the 

. «' 

covariance between the multipole spectra. Gov . In lin- 
ear regime, the covariance neglecting the shot noise con- 
tribution is determined by the galaxy power spectrum in 
real space and parameter /3 = //6, and the off-diagonal 

02 20 

component Gov = Gov is roughly proportional to /3. 
Thus, as increasing the clustering bias h while keeping the 

02 

growth-rate parameter, the covariance Gov becomes 
smaller, and the monopole and quadrupole power spec- 
tra become statistically independent. To see this more 
explicitly, we define 



Gov 



-0,2 



0,0^ 2,2„ • 

[Gov Gov ]i/2 



(20) 



In Fig. m taking account of the shot noise contribution, 
the quantity Tcov is plotted against the parameter /3. 

i-n 

Here, we used the linear theory to calculate Gov . Fig.|3] 
implies that in our fiducial setup with / = 0.858, Tcov 



becomes < 0.2 for the bias 6 = 4. Since the smaller val- 
ues of /3 also suppress the Kaiser effect in the covariances 

00 22 

Gov and Gov , the constraints from the monopole and 
quadrupole spectra is relatively improved. 

The result suggests that even the partial information 
with monopole and quadrupole spectra still provides a 
fruitful constraint on Da-, H and /, depending on the 
survey setup. In this respect, a benefit to use these power 
spectra should be further explored. As a next step, we 
will discuss the robustness of the parameter constraints 
against the systematic biases. 



C. Impact of systematic biases 

Among various envisaged systematics that affect the 
parameter constraints, the incorrect assumption for the 
theoretical template of power spectra may seriously lead 
to a bias in the best-fit parameters. There are several 
routes to produce an incorrect theoretical template; in- 
correct model of redshift distortion and/or non-linear 
gravitational evolution, wrong prior information for cos- 
mological parameters, and improper parametrization for 
galaxy bias. In this subsection, we specifically examine 
the first and second cases. We first discuss the incorrect 
model of redshift distortion, and quantify the size of the 
systematic bias in the best-fit parameter. The influence 
of the wrong prior information will be discussed in next 
subsection. 



1. Systematic biases from a wrong model of redshift 
distortion 

Let us first discuss the impact of incorrect model of 
redshift distortion on the parameter estimation. To be 
precise, we consider the small discrepancy in the theoret- 
ical template for redshift-space power spectrum (jlSp , and 
estimate the systematic biases from Eq. pTj) . Fig.[S]shows 
the systematic biases caused by the incorrect model tem- 
plate neglecting the A and B terms. We plot the re- 
sults by varying the model parameters, A:i„ax (top- left), 
Ug (top-right), b (bottom- left), and ctv (bottom-right), 
around the fiducial values. In each panel, the first three 
panels from the top plot the deviation of the best-fit value 
from the fiducial one, Sf, SDa, and 6H, normalized by 
their fiducial values. On the other hand, the lowest panel 
shows the Figure-of-Bias (FoB), which represents the sta- 
tistical significance of systematic biases relative to the 
statistical errors, defined by [47ll48j: 



1/2 



FoB = ij2 SO^F', 



(21) 



Note that the matrix F'^ is the same inverse of the sub- 



matrix i^j as defined in Eq. (|T9|) . but with the Fisher 
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FIG. 5: Systematic biases for best-fit values of parameters /, Da and H and Figure-of-Bias as function of fcmax (top- left), 
ng (top-riglit), h (bottom-left), and CTv (bottom-right). These are the estimates adopting the 'incorrect' model of redshift- 
space power spectrum, in which we ignore the small correction terms, A and B. In bottom panels of each figure, thick and 
thin lines respectively show the FoB in three and two-dimensions, i.e., {Da, H , f) and {Da,H). The dotted lines indicates 
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unchanged irrespective of the survey volume Va, while the FoB given here represents the specific results with the survey volume 
Vs = 4/i~^Gpc''. The fiducial values of the model parameters used in the calculation are the same as in Fig. [3] (indicated by 
vertical dotted lines), except the parameter along the horizontal axis. 



matrix obtained from the incorrect template. With the 
above definition, the FoB squared simply reflects the Ax^ 
for the true values of the parameters relative to the bi- 
ased estimate of the best-fit values [i^. Thus, in the 
cases with three parameters, if the FoB exceeds 1.88 (in- 



dicated by the red, thick dotted lines), the true values 
of the parameters would go outside the l-cr (68%C.L.) 
error ellipsoid of the biased confidence region. Notice 
that the shift of best-fit parameters remains unchanged 
irrespective of the survey volume Vg, while the FoB is 
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proportional toVs ■ 

Fig. [5] shows that the biases in the distance informa- 
tion, 8Da and Jff, are basicaUy small and reach 1 — 2% 
at most, but the bias in the growth-rate parameter, (5/, 
is rather large. Hence, the behaviors of the FoBs indi- 
cated by the thick lines are mostly dominated by the 
error and bias in the growth-rate parameter. As a result, 
for some ranges of parameters, the expected FoB using a 
full-shape information (black solid, labeled as 'full 2D') 
tends to exceed the critical value, 1.88. This is true even 
if we marginalize over / and just focus on the distance in- 
formation Da and -ff , depicted as thin lines in the lowest 
panels (labeled as 'FoB2d')- Note that in the case of two 
parameters, the true values of and H are ruled out 
at 1-cr level if FoB exceeds 1.52 (red, thin dotted lines). 

On the other hand, if we use the information ob- 
tained only from the monopole and quadrupole spectra 
(magenta, dot-dashed lines), the systematic biases are 
significantly reduced, and the resultant FoBs are well 
within the critical values except for unrealistic case with 
a large Uv or anti-bias h < \. If we are just inter- 
ested in Da and H marginalized over /, the FoB be- 
comes substantially smaller, and would be far below the 
critical value 1.52, even for a large galaxy survey with 
Vs < 40/i~'^Mpc^. Therefore even the partial informa- 
tion from the monopole and quadrupole spectra is help- 
ful and rather robust against the systematic biases than 
the full 2D information. Although the figure-of-merit for 
the constraints on Da-, H and / would be degraded, the 
reduction of FoM is at most factor of ~ 0.6, which can be 
improved to ^ 0.8 for highly biased objects (see Fig. 

Finally, there are several interesting points to be noted. 
One is the oscillatory behavior of the systematic biases 
and FoB shown in the top-left panel. This is originated 
from the acoustic structure of the power spectrum, and 
the result suggests that the bias in the growth-rate pa- 
rameter 5f is sensitively affected by the BAO measure- 
ment. Another noticeable feature is a suppression of the 
FoB in the case of three parameters using the monopole 
and quadrupole spectra, which appears at a larger value 
of the galaxy bias b (thick, dot-dashed line in bottom- 
left panel). This is presumably due to the multiple ef- 
fects that as increasing the clustering bias, the systematic 
bias for the growth-rate parameter tends to be slightly re- 
duced, while the constraint on the growth-rate parameter 
becomes gradually weaker. There also appears a similar 
trend in the case using a full 2D spectrum, but the sup- 
pression is rather small and FoB never falls below the 
critical value, 1.88. This is because the biased estimate 
of growth-rate parameter, 5f, significantly deviates from 
the fiducial value, as opposed to the case using monopole 
and quadrupole spectra. 



2. Systematic biases from incorrect prior information 

So far, we have assumed that the underlying cosmolog- 
ical parameters necessary to compute the redshift-space 



power spectrum are whole known a priori from the CMB 
observations such as PLANCK. However, even the preci- 
sion CMB measurement produces some uncertainties in 
the cosmological parameters due to the parameter degen- 
eracy. This may give an important source for the incor- 
rect theoretical template for redshift-space power spec- 
trum, and leads to a biased estimate of Da, H, and /. 

Fig. ini quantifies the size of systematic biases and FoB 
arising from the incorrect assumptions for cosmological 
parameters. Here, we especially focus on the parameters 
As, ilm, and h fixing fim^^ constant, and plot the sen- 
sitivity of the systematic biases to the variation of those 
parameters. Note that in computing the power spectrum, 
we strictly assume the flat cosmological model and the 
model of redshift distortion (|15p as a fiducial model of 
power spectrum template. 

Compared to the results in Sec. IIV C 1[ the systematic 
bias in the growth-rate parameter is relatively small, and 
the significance of the biases in the acoustic-scale infor- 
mation conversely increases. That is, the best-fit values 
of the parameters Da and H is rather sensitive to the 
precision of the prior information in the power spectrum 
template. A noticeable point is that this is true irre- 
spective of the choice of the template power spectra used 
in the parameter estimation (i.e., full 2D spectrum or 
combination of Pq and ^2)- As a result, a percent-level 
precision is generally required for the prior information 
of cosmological parameters, except for the scalar spectral 
amplitude. As- Through the non- linear clustering and/or 
redshift distortion, a small change in As alters the power 
spectrum shape, and it can potentially affect the acous- 
tic scale and the clustering anisotropics. However, at 
z = 1, the non-linear effects on the scales of our interest, 
k < 0.2/iMpc^^, is rather mild, and the resultant impact 
on the acoustic-scale measurement is extremely small. 
Hence, for a typical survey volume of stage Ill-class sur- 
vey with Vs ~ 4/i~'^Gpc^, no appreciable systematic bias 
might be produced from the incorrect prior assumption 
on As- 



V. SUMMARY 

In this paper, we have studied the cosmological con- 
straints from the anisotropic BAOs based on the mul- 
tipole expansion of redshift-space power spectrum. We 
have derived the several formulae for the Fisher analysis 
using the multipole power spectra; Eqs. ([7]) and © for 
the Fisher matrix, and Eqs. (ITT]) and for the esti- 
mation of systematic biases. We then consider the hypo- 
thetical galaxy survey of Vs — Ah^^Gpc^ and z = 1, and 
discuss the potential power of the lower multipole spec- 
tra on the cosmological constraints, particularly focusing 
on the parameters Da, H and /. 

Compared to the analysis with full 2D power spectrum, 
a partial information from the monopole and quadrupole 
power spectra generally degrades the constraints on Da, 
H, and /. Typically, the constraint is degraded by a fac- 
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FIG. 6: Systematic biases for the best-fit values of the parameters /, Da and H , and FoB for these three parameters (from top 
to bottom), adopting the incorrect prior information for cosmological parameters in computing the template power spectrum; 
X — As, and h {Q.mh'^- fixed). The results are plotted against the fractional difference between the correct and incorrect 
values of each cosmological parameters, AX/Xnd- Solid and dashed lines represent the results from a full 2D power spectrum 
and partial information with monopole and quadrupole spectra, respectively. Note that in bottom panel, the horizontal dotted 
lines indicates the 1-a significance of the deviation relative to the statistical error. 



tor of ~ 1.3 for each parameter. The interesting finding is 
that adding the information from hexadecapole spectra 
(Pi) to that from the monopole and quadrupole spec- 
tra greatly improves the constraints, and the resultant 
constraints would become almost comparable to those 
expected from the full 2D power spectrum (see Fig. [3]) . 
Note also that the situation would be relatively improved 
depending on the properties of galaxy samples, and for 
highly biased galaxy samples with 6 ~ 4, the total power 
of the constraints defined by the Figure-of-Merit [FoM, 
Eq. can reach ~ 80% of the one expected from the 
full 2D power spectrum. 

We have also investigated the impacts of systematic 
biases on the best-fit values of Da, H and /. The incor- 
rect model of redshift distortion tends to produce a large 
systematic bias in the growth-rate parameter, and the 
size of biases would be rather significant for the analysis 
with full 2D spectrum. An interesting suggestion is that 
the situation would be greatly relaxed if we only use the 
combination of monopole and quadrupole spectra, and 
the estimated value of Figure-of-Bias defined by Eq. ([2T|) 
is mostly below the critical value for stage-Ill class sur- 
veys (Fig. [S]). In this respect, the analysis with partial 
information from monopole and quadrupole may be still 



helpful in cross-checking the results derived from the full 
2D power spectrum. On the other hand, wrong prior as- 
sumption of cosmological parameters in computing the 
template power spectrum severely affects the acoustic- 
scale determination, and a percent-level precision is re- 
quired for the prior information in order to avoid a large 
systematic biases on Da and H (Fig. This is true 
irrespective of the choice of template power spectra used 
in the analysis. 

Finally, we note that the assumptions and situations 
considered in the paper are somewhat optimistic or too 
simplistic, and a more careful study is needed for a 
quantitative parameter forecast. One critical aspect is 
the modeling of the galaxy power spectrum. In real- 
ity, the assumption of linear and deterministic galaxy 
biasing is idealistic, and the scale-dependence or non- 
linearity/stochasticity of the galaxy biasing should be 
consistently incorporated into the theoretical template of 
redshift-space power spectrum. Although this is tiny ef- 
fect for the scale of our interest, the distance information. 
Da and H, is rather sensitive to a slight modification of 
the acoustic structure in the power spectrum, and results 
in this paper might be somehow changed. A more elabo- 
rate modeling for power spectrum is thus quite essential. 



12 



Acknowledgments 

AT is supported by a Grant-in-Aid for Scientific Re- 
search from the Japan Society for the Promotion of Sci- 
ence (JSPS) (No. 21740168). TN is supported by JSPS. 
SS is supported by JSPS through research fellowships 
and Excellent Young Researchers Overseas Visit Pro- 



gram (No. 21-00784). This work was supported in part 
by Grant-in-Aid for Scientific Research on Priority Areas 
No. 467 "Probing the Dark Energy through an Extremely 
Wide and Deep Survey with Subaru Telescope", JSPS 
Core-to-Core Program "International Research Network 
for Dark Energy" , and World Premier International Re- 
search Center Initiative (WPI Initiative), MEXT, Japan. 



D. J. Eisenstein et al. (SDSS), Astrophys. J. 633, 560 [29 
(2005), astro-ph/0501171. 

W. J. Percival et al., Mon. Not. Roy. Astron. Soc. 381, [30 
1053 (2007), 0705.3323. 

W. J. Percival et al., Mon. Not. Roy. Astron. Soc. 401, 

2148 (2010), 0907.1660. [31 

C. Alcock and B. Paczynski, Nature 281, 358359 (1979). 

H.-J. Seo and D. J. Eisenstein, Astrophys. J. 598, 720 [32 

(2003), astro-ph/0307460. [33 

C. Blake and K. Glazebrook, Astrophys. J. 594, 665 [34 
(2003), astro-ph/0301632. [35 
M. Shoji, D. Jeong, and E. Komatsu, Astrophys. J. 693, 

1404 (2009), 0805.4238. [36 
N. Padmanabhan and M. J. White, 1, Phys. Rev. DTT, 
123540 (2008), 0804.0799. [37 

E. V. Linder, Astropart. Phys. 29, 336 (2008), 0709.1113. 
L. Guzzo et al., Nature 451, 541 (2008), 0802.1944. [38 
K. Yamamoto, T. Sato, and G. Huetsi, Prog. Theor. 
Phys. 120, 609 (2008), 0805.4789. [39 
Y.-S. Song and W. J. Percival, JCAP 0910, 004 (2009), 
0807.0810. [40 

D. Schlegel, M. White, and D. Eisenstein (with input 

from the SDSS-III) (2009), 0902.4680. [41 

G. J. HiU et al. (2008), 0806.0183. [42 
J. P. Beaulieu et al. (2010), 1001.3349. 
N. Gehrels (2010), 1008.4936. 

H. -J. Seo and D. J. Eisenstein, Astrophys. J. 665, 14 [43 
(2007), astro-ph/0701079. 

M. White, Y.-S. Song, and W. J. Percival, Mon. Not. [44 
Roy. Astron. Soc. 397, 1348 (2008), 0810.1518. [45 
T. Okumura et al., Astrophys. J. 676, 889 (2008), 
0711.3640. ' [46 

R. Takahashi et al., Astrophys. J. 700, 479 (2009), [47 
0902.0371. 

A. Cabre and E. Gaztanaga, Mon. Not. Roy. Astron. Soc. [ 
393, 1183 (2009), 0807.2460. ' [49 

N. Kaiser, Mon. Not. Roy. Astron. Soc. 227, 1 (1987). 
A. J. S. Hamihon, Astrophys. J. 385, L5 (1992). 
A. J. S. Hamihon (1997), astro-ph/9708102. 

D. Tocchini-Valentini, M. Barnard, C. L. Bennett, and [50 
A. S. Szalay (2011), 1101.2608. 

E. Komatsu et al. (WMAP), Astrophys. J. Suppl. 180, 

330 (2009), 0803.0547. [51 
K. Yamamoto, Astrophys. J. 595, 577 (2003), astro- 
ph/0208139. 

M. Tegmark, Phys. Rev. Lett. 79, 3806 (1997), astro- 
ph/9706198. 



H. A. Feldman, N. Kaiser, and J. A. Peacock, Astrophys. 

J. 426, 23 (1994), astro-ph/9304022. 

K. Yamamoto, M. Nakamichi, A. Kamino, B. A. Bassett, 

and H. Nishioka, Publ. Astron. Soc. Jap. 58, 93 (2006), 

astro-ph/0505115. 

S. Saito, M. Takada, and A. Taruya, Phys. Rev. D80, 
083528 (2009), 0907.2922. 

A. Taruya, T. Nishimichi, and S. Saito (2010), 1006.0699. 

R. Takahashi et al. (2009), 0912.1381. 

M. C. Neyrinck and I. Szapudi (2007), 0710.3586. 

A. Taruya, T. Nishimichi, S. Saito, and T. Hiramatsu, 

Phys. Rev. D80, 123503 (2009), 0906.0507. 

D. J. Eisenstein and W. Hu, Astrophys. J. 496, 605 

(1998), astro-ph/9709112. 

R. Scoccimarro, Phys. Rev. D70, 083007 (2004), astro- 
ph/0407214. 

M. Crocce and R. Scoccimarro, Phys. Rev. D77, 023533 
(2008), 0704.2783. 

T. Okumura and Y. P. Jing, Astrophys. J. 726, 5 (2011), 
1004.3548. 

D. Jeong and E. Komatsu, Astrophys. J. 691, 569 (2009), 
0805.2632. 

S. Saito, M. Takada, and A. Taruya (2010), 1006.4845. 
W. E. Ballinger, J. A. Peacock, and A. F. Heavens, 
Mon. Not. Roy. Astron. Soc. 282, 877 (1996), astro- 
ph/9605017. 

H. Magira, Y. P. Jing, and Y. Suto, Astrophys. J. 528, 

30 (2000), astro-ph/9907438. 

A. J. Albrecht et al. (2006), astro-ph/0609591. 

A. Taruya and T. Hiramatsu, Astrophys. J. 674, 617 

(2008), 0708.1367. 

(2006), astro-ph/0604069. 

S. Joudaki, A. Cooray, and D. E. Holz, Phys. Rev. D 80, 
023003 (2009), 0904.4697. 

C. Shapiro, Astrophys. J. 696, 775 (2009), 0812.0769. 
Throughout the paper, we work with the distant-observer 
approximation, and neglect the angular dependence of 
the line-of-sight direction, relevant for the high-redshift 
galaxy surveys. 

Here, we use the standard notation for the multipole ex- 
pansion of redshift power spectra given by ([1]), which 
differs from the definition of Ref. [SQl 
The definition of velocity dispersion CTv adopted in this 
paper differs from the one commonly used in the lit- 
erature by a factor of /, but coincides with those in 
Refs. 



